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Abstract 

By Monte Carlo simulations of a variant of the bond-fluctuation model without topological 
constraints we examine the center-of-mass (COM) dynamics of polymer melts in d = 3 dimen- 
sions. Our analysis focuses on the COM displacement correlation function C^{t) ~ dfh^{t)/2, 
measuring the curvature of the COM mean-square displacement h-t<}{t). We demonstrate that 
CN(t) ~ -{R^/T^f{p*/p) f{x = t/TN) with N being the chain length (16 < < 8192), 
i?N ~ N^^'^ the typical chain size, Tn ~ A^^ the longest chain relaxation time, p the monomer 
density, p* ^ N/R^ the self-density and f{x) a universal function decaying asymptotically as 
f{x) ~ x~'^ with oj = {d + 2) X a where q = 1/4 for x <^ 1 and a = 1/2 for x ^ 1. We argue that 
the algebraic decay NC-^{t) ~ — for t ^ Tn results from an interplay of chain connectivity 
and melt incompressibility giving rise to the correlated motion of chains and subchains. 
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Electronic address: joachim.wittmer@ics-cnrs.unistra.fr 





u(0) 




insity 












M 

t=ol 





FIG. 1: The displacement auto-correlation function C{t) of dense colloids without momentum 
conservation is known to reveal a negative algebraic tail due to long-range dynamical interactions 
caused by the weak compressibility of the solution [2]: If a particle (filled sphere) is displaced 
at time t = by a distance m(0) along the x-axis, this creates a density dipole (bold peaks) 
decaying in time by the cooperative diffusion of the density field (dashed line). The gradient of 
the density field generates a force pulling the particle (open sphere) back to its original position. 
We argue that a related mechanism is relevant for describing the correlations of the chain and 
subchain center-of-mass displacements in dense polymer melts without topological constraints and 
momentum conservation. 

I. INTRODUCTION 



Prelude: Overdamped colloidal suspensions. Dense, essentially incompressible simple 
liquids with conserved momentum are known to exhibit long-range correlations of the particle 
displacement field [l-sl. As first shown in the seminal molecular dynamics (MD) simulations 
by Alder and Wainwright js], the coupling of displacement and momentum fields manifests 
itself by an algebraic decay of the velocity correlation function (VCF), 



C{t) = {v{t) ■ v{0)) ~ +l/e(t)'' ~ 



ad 



(1) 



with v{t) being the particle velocity at time t, d the spatial dimension and ^{t) ~ t° the 
typical particle displacement with exponent a = 1/2 j4|. Interestingly, even if the momentum 
conservation is dropped, as justified for overdamped dense colloidal suspensions Q], scale- 
free albeit much weaker correlations of the displacement field are to be expected due to the 
incompressibility constraint As illustrated in Fig. [H the motion of a tagged colloid is 



2 



coupled to the collective density dipole field 

, , , 1 u(0)-r 



e~ir/mf for t > 0, 



(2) 



created by the colloid's own displacement n(0) at t = 0. After averaging over the typical 
displacements of the test particle and assuming a Cahn-Hilliard response proportional to the 
gradient V/i(r) of the chemical potential fi{r) of the density field, this leads to a negative 
algebraic long-time decay of the VCF 



d+2 



-Ijr with w = (c? + 2) a, 



(3) 



e.g., Ci; = 5/2in(i = 3 dimensions. This phenomenological scaling picture agrees with mode- 
coupling calculations 0,0, Q| and has been confirmed computationally by means of Lattice- 
Boltzmann simulations Lfll, MD simulations ISll and even Monte Carlo (MC) simulations with 
local moves (as discussed below) [9[. 

Deviations from Flory's ideality hypothesis. In this study we explore the dynamics of 
another complex fluid for which momentum conservation is generally believed to be irrel- 
evant melts of long and flexible hornopolymers of length N in d = 3 dimensions. 



Following Flory's "ideality hypothesis" 10|, lUj , one expects these chains to obey a Gaussian 
statistics with a typical chain size i?N ~ hN^^"^ where h stands for the effective bond length 
of asymptotically long chains 10|. Recently, this cornerstone of polymer physics has been 
challenged both theoretically and numerically for three-dimensional (3D) melts isj-jj], for 
effectively two-dimensional (2D) ultrathin films 15|, ll6| and one- dimensional (ID) thin cap- 



illaries 



18| . The physical idea behind the predicted long-range correlations is related 
to the "segmental correlation hole", p*{s) ~ s/R'^, of a subchain of arc- length s of typical 



size Rs 
U*{s) ^ 
together 
U*(s) ~ 



. Due to the overall incompressibility of the melt this sets an entropic penalty 
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with p being the total monomer density) against bringing two subchains 



151] . In d = 3 dimensions, the segmental correlation hole effect is weak. 



decay 



12 



l/i/i, and a perturbation calculation can be performed 12|, |l3|. The detailed 

19| an algebraic 



calculation yields, e.g., for the intrachain angular correlations function Pg 



p. 



for 5( < s < A^, 



(4) 



^vr^ b^p s'^/"^ 

at variance to Flory's hypothesis with g being the subchain length spanning the screening 
length of the density fluctuations ^]. Since Pg ~ "^fPg, the angular correlation function 



provides a direct measure of the curvature of the mean-squared sub chain size i?g and Eq. 
imphes that ^ U*{s) ~ lU/s 



compressibihty kt of the solution 14 
total monomer structure factor 
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13| . Interestingly, g is related to the isothermal 



20| and can thus be determined directly from the 



1 



N-^oo 



g = T ktp, 



(5) 



ra,m=l 



r„ being the position of monomer n, q the conjugated wavevector, V = L'^ the volume of 
the system and T the temperature. (Boltzmann's constant is set to unity throughout 
this paper.) Due to its definition, g is often called "dimensionless compressibility" jisl . 
Remarkably, Eq. does no depend explicitly on g. This reflects the fact that the deviations 
arise due to the incompressibility of "blobs" 11 1 on scales corresponding to s ^ g Q. 



Aim of this study. Naturally, Eq. (jl]) and related findings beg the question of whether a 
similar interplay between the connectivity of the chains and the incompressibility of the melt 
may cause measurable scale-free and gf- independent dynamical correlations between chains 
and between subchains. To avoid additional physics and to simplify the problem we focus on 
polymer melts where hydrodynamic 2l| and topological constraints may be considered to be 



negligi 
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e, as in the pioneering work by Paul et al. 



22 



24 



23l . or are deliberately switched off 

n 

25| . Deviations from the expected Rouse- type dynamics [10[ have indeed been 



reported for such systems in various numerical 22 
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29 
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and experimental studies 
3l|. Characterizing, e.g., the motion of the chain center-of-mass (COM) r-^{t) by its 



mean-square-displacement (MSD) h-^{t), it was found that 

h^{t) = {{r^{t) - rN(0))2> ~ t^ for t < Tn 



(6) 



with Tn ~ A^^ being the longest (Rouse) relaxation time and f3 ^ 0.8 an empirical exponent. 
This is of course at variance to the key assumption of the Rouse model that the random 
forces acting on the monomers (and thus on the chain) are uncorrelated, which implies an 
exponent {3 = 1 for all times [lOj . In this paper we attempt to clarify this problem by means 
of MC simulations of a variant of the bond- fluctuation model (BFM) {32] using local hopping 
moves which do not conserve topology 3] and assuming a finite monomer overlap penalty 
e [l^. In analogy to our work on the angular correlation function Ps ~ d^R^ 12], our 
analysis focuses on the COM displacement correlation function Cj<f(t) ~ d^hj<f(t), measuring 
the curvature of /iN(^) M- This allows us to probe directly the colored forces from the 



molecular bath [2| acting on a tagged reference chain displaced at t = and, in turn, to 
elaborate the scaling theory sketched below. 

Key results. We demonstrate numerically that the VCF C^{t) does indeed not vanish 



lOj. Instead, 



as it would if all the random forces acting on the chains were uncorrelated 
CT<f(t) is found to scale as 

c.w.-(^!i)^^!W/,./r.) (7) 

with /(x) being a universal scaling function. Note that the postulated Eq. ([7]) does not 
depend explicitly on the compressibility of the solution. The squared characteristic "veloc- 
ity" (/^n/^n)^ arises for dimensional reasons. The prefactor p*{N)/p is motivated by the 
correlation hole penalty, i.e. the incompressibility constraint which ultimately couples the 



displacements of (sub)chains [33|. As one expects from Eq. ([3]), the scaling function decays 
as f{x) ~ l/x"^ with an exponent u = {d + 2)/2 for x ^ 1. We show that this long-time 



behavior is preceded by a much weaker algebraic decay with an exponent 34 1 



w = (rf + 2)a = 5/4 for X < 1 (8) 

due to the much slower relaxation, a = 1/4, of the collective dipole field of subchains which 
was generated by the initial displacement of a tagged subchain at t = 0. The gradient of the 
chemical potential V/x(r) pulling the reference subchain back to its original position is of 
course not only due to the density fluctuation of the subchain density field but also to the 
tensional forces along the chains caused by the displacement. (Subchain density fluctuations 
and tensions are coupled and associated both to a free energy fluctuation of order k-^Ts^N'^.) 
Since C^(t) ~ d^h-^it), it follows from Eq. (IHl) that for sufficiently short times, such that the 
white forces acting on the chains are negligible, we expect to find 

/3 = 2- w = (6- rf)/4 = 3/4, (9) 



which is rather similar to the fit, Eq. ([8]), suggested in the literature. See Refs. [35| and [36| 



for two close 
approach 



y related theoretical studies. Interestingly, Schweizer's mode-coupling theory 
35l | is consistent with Eq. We stress that the presented numerical study is 
necessarily incomplete since important dynamical correlations are expected to arise in more 
reahstic models due to topological constraints and, even more importantly, due to not fully 
screened hydrodynamic interactions which have recently been shown to matter [21,] . 



Outline. The paper is organized as follows. In Sec. [IT] the numerical algorithm is intro- 
duced and some technical details are discussed. Our computational results are presented in 
Sec. Illllwhere we focus on essentially incompressible melts (Sec. IIII Bp but also discuss effects 
of finite excluded volume (Sec. IIII C\i . Our results are summarized in Sec. IIV A[ The paper 
concludes in Sec. IIVBI with a comment on what we would expect if topology conservation 
is switched on again. 



II. SOME ALGORITHMIC DETAILS 
A. The classical bond-fluctuation model 



The classical BFM is an efficient lattice MC algorithm for coarse-grained polymer chains 
where each monomer occupies exclusively a unit cell of 2^^ lattice sites on a rf- dimensional 
simple cubic lattice 22|, |23|, |32[. (The fraction (p of occupied lattice sites is thus (p = 2'^p.) 



The BFM was proposed in 1988 by Carmesin and Kremer [32l] as an alternative to single-site 
self-avoiding walk models, which retains the computational efficiency of the lattice without 
being plagued by ergodicity problems. The key idea is to increase the size of the monomers 
and the number of bond vectors to allow a better representation of the continuous-space 
behavior of real polymer melts. A widely used choice of bond vectors for the 3D variant 
of the BFM is given, e.g., by all the permutations and sign combinations of the six vectors 



22,123 



26 



27| 



(2, 0, 0), (2, 1, 0), (2, 1, 1), (2, 2, 1), (3, 0, 0), (3, 1, 0). 



(10) 



If only local MC moves of the monomers to the six nearest neighbor sites are performed — 
called "L06" moves [[ISi] — this vector set ensures automatically that polymer chains cannot 
cross. (These L06- moves are represented in Fig. MJo) by the filled circles.) Consequently, 
several authors report a reptation-type dynamics for chain lengths above A^ ~ 200 at a 



"melt" volume fraction 



ip = 0.5 [22 



23 



26 



27|. 



B. BFM with topology violating local moves 

For consistency with these studies we keep Eq. ( ITOj) . although the non-crossing constraint 
is irrelevant for us, and use the same volume fraction (p = 0.5. With respect to the classical 



6 




108 allowed bonds 



FIG. 2: The BFM is a lattice MC algorithm for coarse-grained polymer chains where monomers 
are represented by cubes on a simple cubic lattice connected by a set of allowed bond vectors 



(a) We use here a variant of the BFM where a finite energy e has to 



overlap and a corresponding fraction for a partial monomer overlap 



3e paid if two cubes totally 



]J]. (b) Using local MC jump 



attempts to the next (filled circles) and next-nearest (open circles) neighbors we investigate the 
influence of the incompressibility constraint on the dynamics of polymer melts without topological 
constraints, (c) In Sec. IIII Al we will consider beads {N = 1) confined between hard walls [igI. L J| . 



variant of the BFM in c? = 3 our algorithm differs in two important points: 

(i) As described in Ref. jl4l], we use a finite excluded volume penalty e which has to be 
paid if two monomers fully overlap. The overlap of two cube corners is sketched in Fig. [2](a). 
Temperature is arbitrarily set to unity. Technically, the finite monomer interaction pen alty 
is implemented using a Potts spin representation of the discrete local density field p(r) 14 1. 
The monomer overlap penalty parameter e is in fact a Laplace multiplier controlling the 
fluctuations of spins and imposing thus the incompressibility of the melt. By reducing the 
fluctuations of the spins the Laplace multiplier causes thus the effective entropic forces which 
imply the dynamical correlations discussed in Sec. IIIII This modification allows us to check 
at one constant volume fraction (0 = 0.5) that the dynamical correlations do not depend 
explicitly on the monomer excluded volume as stated in Eq. (JTj). 

(a) We use local hopping moves to the 26 next and next-nearest lattice sites (so-called 
"L26" moves). Due to the larger jumps, chains cross and the dynamics is of Rouse-type: 
the chain self-diffusion coefficient D-^, e.g., scales as ~ for all chain lengths even 
if monomer overlap is disallowed {e = oo) [3]. 
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0.0 


oo 


2.718 


2.72 


0.2109 


0.032 


0.065 


0.01 


209 


2.718 


2.80 


0.2109 


0.030 


0.062 


0.1 


22 


2.719 


2.92 


0.2067 


0.024 


0.058 


1 


2.4 


2.721 


3.13 


0.1796 


0.015 


0.040 


10 


0.32 


2.670 


3.24 


8.8E-02 


0.003 


0.009 


100 


0.25 


2.636 


3.24 


6.9E-02 


0.0010 


0.003 



TABLE I: Various properties for BFM polymers melts of arbitrarily large chain length A'^ at volume 
fraction (p = 8p = 0.5: the dimensionless compressibility g, the root-mean-square bond length I, 
the effective bond length b, the acceptance rate A, the local monomer mobility W and the self- 
diffusion coefhcient Dn- The dynamical data refer to L26-moves to the nearest and next-nearest 
lattice sites. 

Although the use of L26-moves and finite monomer interactions does speed up the re- 
laxation dynamics, it remains obviously impossible to equilibrate dense polymer solutions 
with chain lengths up to A^ = 8192 just using local hopping moves. Taken advantage of our 
previous studies on static properties Q, Q, we have used configurations equilibrated using 
a mix of global slithering snake and double bridging moves together with local L26-moves. 
We use periodic simulation boxes of linear dimension L = 256. Thus at = 8p = 0.5 these 
systems contain NM = pL^ = 2^° ^ 10^ monomers and even for chains with A^ = 8192 we 
still have M = 128 chains. 

Using these well characterized configurations we have computed time series over 10'' or 
10® Monte Carlo Steps (MCS) using local L26-moves. We stress that the aim is not to cover 
necessarily the full relaxation time Tn ~ A^^ for our larger chains, but to precisely describe 
the dynamical correlations which are the most pronounced at short times. As main working 
point we take e = 10 for which the chain length has been scanned up to A^ = 8192. For 
other penalties e we present results obtained using A^ = 1024. Some relevant properties are 
summarized in Table [T] and are represented in Fig. [3] and Fig. HI 
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FIG. 3: Summary of static properties for BFM chains at volume fraction (p = 8p = 0.5: (a) 
Root mean-square bond length l{e) (triangles) and effective bond length b{e) (squares) in d = 3 
dimensions, (b) The dimensionless compressibility g{e) for asymptotically long chains levels off for 
large e where g{e S> 10) ~ 0.25 (bold line), (c) Pg for several e compared to the exponent —3/2 
(bold line) expected for s » (7 according to Eq. (jH [ij] . 

C. Reminder of static properties 



Panel (a) of Fig. [3] shows the mean-square bond length 1(e) and the effective bond length 
b^e) for asymptotically long chains as determined in Ref. 1^. Obviously, b{e) — )■ l{e = 
0) = 2.718 in the small-e limit. b{e) then increases in the intermediate e-window before it 
levels off at our main working point e = 10. As one may expect, systems with e = 100 
cannot be distinguished from systems computed using the classical BFM without monomer 
overlap [e = 00). Panel (b) presents the dimensionless compressibility g{e) obtained from 
total static structure factor, Eq. ([5]). The dashed line corresponds to the asymptotics for 
weak interactions, l/g{s) ~ e [14]. The compressibility levels off for large penalties where 
g{6 ^ 10) — 0.25 (bold line). The angular correlation function Pg 19[ is presented for 
several overlap penalties e and one chain length = 1024 in panel (c). The theoretical 
prediction Eq. (|4]) is nicely confirmed for s ^ g{e) by the power-law asymptote indicate by 
the bold line. 
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FIG. 4: Acceptance rate A, effective local mobility W and self-diffusion coefficient ND^ obtained 
using L26-moves at volume fraction (f) = 0.5 vs. the monomer overlap penalty e. All dynamical 
properties decrease monotonously with the interaction penalty but become essentially constant 
above our main working point at e = 10. The dashed line corresponds to the diffusion coefficient 
according to the Rouse prediction, Eq. (fT3]l . assuming the indicated mobilities W (squares). 



III. COMPUTATIONAL RESULTS 



A. Diffusion of dense BFM beads 



Since in MC simulations there is no "monomer mass", no (conserved or non-conserved) 
"monomer momentum" and not even an instantaneous velocity, it might at first sight appear 
surprising that a well-posed "velocity correlation function" (VCF) can be defined and mea- 
sured. To illustrate that this is indeed the case is the first purpose of this subsection. The 
second is to verify that the negative analytic decay of the VCF expected for overdamped 
colloids, Eq. ([3]), is also of relevance for dense BFM beads (A^ = 1) diffusing through configu- 
ration space by means of local hopping moves on the lattice. The systems presented in Fig. [5] 
correspond to three different effective dimensions d. The effectively ID systems {d = 
have been obtained by confining the beads to a thin capillary of square cross-section, the 2D 
systems {d = 2+) by confining the beads to a thin slit as shown in Fig. [2]^c). The distance 
H = 4 between parallel walls allows the free crossing of the beads. (For the ID case this is 
crucial, of course.) The data has been obtained for beads without monomer overlap {e = oo) 
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t [MCS] 

FIG. 5: Diffusion of BFM beads in d = 1^, d = 2+ and d = 3 dimensions obtained using L26-moves. 
Inset: Altliough the bead diffusion is essentially free, small deviations are visible if h{t)/{2dDt) 
is plotted in log-linear coordinates. Main panel: Collapse of C{t,6t) for d = 1'^ and various time 
increments 6t (open symbols). The "final" function C{t) is obtained by adding the first decade 
of data for each 5t and logarithmic averaging. The cummulants agree nicely with the predicted 
exponent w = (d + 2)/2, Eq. 

by means of L26-moves on a 3D cubic lattice for a volume fraction (p = 0.5. We average 
over the 2^° beads contained in each configuration. 

One standard measure characterizing the monomer displacements is the MSD h{t) = 
{{r(t) — r(0))^) displayed in the inset of Fig. |5] (with r(t) being the particle position at time 
t). As one expects, the displacements become uncorrelated for t ^ 1, i.e. h{t) ~ 2dDt with 
D being the monomer self-diffusion constant: D = 0.0187 for d = 1+, D = 0.0154 for d = 2+, 
D = 0.0382 for d = 3. The typical displacement ^(t) = h^^^(t) thus scales as ^{t) ~ with 
a = 1/2. However, small deviations are clearly visible for short times if h{t)/{2dDt) is 
plotted in log-linear coordinates. The deviations are particulary strong for d = The 
effective random forces acting on the beads are thus not completely white. Obviously, one 
might try to characterize these deviations by fitting various polynomials to the measured 
MSD. Since one needs to substract, however, the huge free diffusion contribution from the 
measured signal to obtain tiny deviations this is a numerical difficult if not impossible route. 

In analogy to the static angular correlation function Ps ~ d^R^ allowing to make mani- 
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fest deviations from the Gaussian chain assumption, it is numerically much better to directly 
compute the second derivative of h{t) with respect to time. How this can be done is illus- 
trated in the main panel of Fig. O We sample equidistant series of configurations at time 
intervals 6t = 1, 10, 100, ... as indicated by the open symbols. Each time series contains 10^ 
configurations. Averaging over all possible pairs of configurations (toi^o + t) we compute 
the displacement correlation function C(t,6t) = {u(tQ + t) ■ uito)) ^^J 5t^ , i.e. a four-point 
correlation function of the monomer trajectories with u[t) = r{t + 6t) — r{t) being the 
monomer displacement vector at time t in a time interval 6t. By construction C{t, 6t) = if 
both displacement vectors are uncorrelated. Note that if one computed times t shorter than 
6t, both displacement vectors would become trivially correlated, since they describe in part 
the same particle trajectory. Hence, C{t = 0,6t) = h{5t)/5t^ and C{t < 6t,6t) > (not 
shown). More importantly, for times t ^ 6t one has 

^ U(t + U)-2m + U(t-M) ^ 1 

2ot^ 2 

as one readily sees by applying finite-difference operators with respect to time to the 
monomer MSD. As can be seen for d = 1+, the 5/!:- dependence drops indeed out for t/6t > 1. 
We thus often avoid the second index 6t and write C{t) for the displacement correlation 
function. Obviously, the statistics deteriorates for large t/6t where fewer configuration pairs 
contribute to the average (taking apart that the signal itself decays). It is for this reason 
that we need a hierarchy of time series of different 6t. Taking for each time series only 
the first decade of data (2 < t/6t < 20), these reduced data sets are pasted together and 
averaged logarithmically. These logarithmic cummulants are given for each dimension and 
compared to the exponents u = 3/2 (thin line), u = 2 (dashed line) and u = 5/2 (bold line) 
predicted by Eq. ([3]) for d = 1, d = 2 and d = 3, respectively. The data agrees over more 
than two orders of magnitude in time with the prediction, especially for d = 1^. If one is 
satisfied with less orders of magnitude it is sufficient to check the exponents using just a time 
window 5t = 1 as can be seen for the open spheres. The superposition of data from time 
series with different 6t is just a numerical trick which reduces the number of configurations 
to be stored and the number of configuration pairs to be computed for a given time t. 
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t [MCS] 

FIG. 6: Various MSD for BFM chains of lengths = 64 and iV = 8192. The open symbols 
refer to the monomer MSD h(t), the filled symbols to the MSD h-^{t) of the COM of chains, the 
crosses to the MSD hs{t) of the COM of subchains of arc- length s = 64 of total chains of length 
N = 8192. The dashed line indicates the monomer MSD expected for Rouse chains for t ^ Tn, 
the thin solid line the free diffusion limit expected for uncorrelated random forces. The bold lines 
represent the effective exponent f3 = 0.8, Eq. The stars show Eq. (|17p . i.e. the time window 
t <C r is described by (3 = 3/4. 

B. Polymer melts without topological constraints 

1. Mean-square displacements 

Having shown that scale-free dynamical correlations exist for dense BFM beads as ex- 
pected for overdamped colloids (Fig. [T]), we turn now our attention to 3D melts of long 
and flexible homopolymers. We focus here on systems with finite monomer overlap penalty 
e = 10 and volume fraction (p = 0.5. Due to the use of the finite overlap penalty and the 
topology non- conserving local L26-moves, the dynamics is indeed of Rouse-type as can be 
seen from the various MSD presented in Fig. [6] for chains of length = 64 (diamonds) and 

= 8192 (spheres). 

The monomer MSD h{t) is indicated by the open symbols. Note that for the large chain 
lengths sampled here it is inessential whether this average is performed over all monomers (as 
in Sec. IIII A|) or only over a few monomers in the (curvilinear) center of chains at n ^ N/2 
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as in Refs. 22|, |23|, |26|. As expected from the Rouse model we obtain for short times 10 1 



h{t) = b^iWty/^N^ for 10 < t < Tn ^ N^W (12) 

as indicated by the dashed hne, i.e. the typical displacement increases as ^{t) ~ with 
a = 1/4. As can be seen for N = 64, the monomers diffuse again freely with a power-law 
slope 1 for times larger than the Rouse time T^. We remind that it was not the aim of 
the present work to sample for our larger chains (A^ > 1000) over the huge times needed 
to make this free diffusion regime accessible. Following Paul et al. 22|, |23|, the short-time 
power law, Eq. fll2p . can be used to determine the effective monomer mobility. We obtain 
W{e = 10) = 0.003 for our main working point. Mobilities for other penalties are listed 
in Tab. HI As may be seen from Fig. HI W{e) decays with increasing excluded volume 
just as the acceptance rate A{e) (spheres) of the hopping attempts, but the decay is even 
more pronounced for W{e), i.e. increasingly more accepted moves do not contribute to the 



effective motion ("cage effect") |23 |. 



The full symbols displayed in Fig. [6] refer to the MSD hj<s{t) of the chain COM defined in 
Eq. (|6]). As one expects for Rouse chains, the amplitude of h^{t) decreases inversely with A^ 
and the diffusion appears to be uncorrelated {h-^{t) ~ t) at least for times t ^ t* ^ 10^ as 
indicated by the vertical arrow. Obviously, the center-of-mass MSD h-^lt) and the monomer 
MSD h(t) merge for times beyond the Rouse time (t ^ Tn). Fortunately, since /in(^) 
becomes linear for t* ^ Tn, it is always possible, even for our largest chain length A^, to 
measure the chain self-diffusion coefficient by plotting h^{t)/Qt in log-linear coordinates. 
For our main working point we thus obtain ND-^{e = 10) ~ 0.009 for all A^. Values for 
other e are again given in Tab. [T] and are represented in Fig. H] (diamonds). That these values 
are consistent with Rouse dynamics can be checked by comparing the measured diffusion 
coefficients with the values obtained using the local mobilities W according to [l^ 

ND^ = l^^ W. (13) 

As shown by the dashed line in Fig. HI Eq. (fT5]l agrees well with the directly measured 
diffusion coefficients. Similarly, it is possible (at least for our shorter chains) to measure the 
longest Rouse relaxation time Tn by an analysis of the Rouse modes and to compare it with 
Tn = AN'^/ir^W [l^. We obtain again a nice agreement between directly and indirectly 
computed relaxation times (not shown). 
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Up to now we have insisted on the fact that our systems are to leading order of Rouse type 
and we have characterized them accordingly. However, deviations from the Rouse picture are 
clearly revealed for short times, especially for h^{t) in agreement with the literature 
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3l|. (That the monomer MSD h{t) also deviates for very short times is due to the lattice 
model, i.e. to a trivial lower cut-off effect associated to the discretization.) In agreement 
with Eq. dnj, the short-time COM motion can be characterized by an effective exponent 
/3 ^ 0.8 (bold lines). Since in our B 



confirms the finding by Shaffer 
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i^M version topological constraints are irrelevant, this 
25l | that the deviations found for short chains using 
the classical BFM algorithm with topological constraints [23] cannot alone be attributed to 
precursor effects to reptational dynamics (as discussed in Sec. lIVBp . 

Before we turn to the more precise numerical characterization of these deviations by 
means of the associated displacement correlation function, let us ask whether the observed 
colored forces acting for short times on the COM of the entire chain are also relevant on 
the scale of subchains of arbitrary arc-length s (1 ^ s < A^). To answer this question we 
compute the MSD displacement hs(t) = ((rs(t) — rs(0))^) associated to the subchain center- 
of-mass rs(t) as shown in Fig. [6] for subchains of length s = 64 in the middle of total chains 
of length = 8192 (crosses). Since for short times the subchain does not "know" that it 
is connected to the rest of the chain, one expects it to behave as a total chain of the same 
length [s = N). This is indeed borne out by our data which are well described by 

sh,{t) ^ Nh^{t) for t <^T,^ s^W (14) 

for all chain length and subchain length s studied. The subchains reveal thus for suffi- 
ciently short times the same colored forces (/3 ~ 0.8) as the total chain as can be clearly seen 
from the example given in Fig. O For larger times, ^ t, the subchain becomes "aware" 
that it is connected to the rest of the chain and gets enslaved by the monomer MSD. We 
thus observe 

h,{t) ^ h{t) ^ h^iWif'^s^N"^ for T,^t<^T^ (15) 
and, obviously, hsit) ^ /in(^) ^ h(t) ^ 6Dn^ for even larger times t T^. 
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2. Locality and relevant exponent a 



Two comments are in order here. First, it should be noticed that Eq. (fT^ expresses 
the fact that the effective forces acting on the N/s subchains of length s in a chain of 
total length N add up independently to the forces acting on the total chain. In this sense 
Eq. ffl^ states that the deviations from the Rouse picture must be local js^]. We will 
explicitly verify this below (Fig. [T0|) . Second, if one chooses following Eq. ffTTl) an arbitrary 
time window 6t to characterize the displacement correlations, this corresponds to dynamical 



blobs containing s ~ VW6t adjacent monomers, which must move together due to chain 
connectivity. Eq. (|T5l) implies now that the dipole field jsSi associated with the COM of these 
s-subchains and created at t = by a tagged s-subchain must decay according to a typical 

1 /2 

displacement C,(t) ^ hs (t) ~ t° with a = 1/4. It is this exponent a which is mentioned 
in the Introduction, Eq. ([8]). Since this exponent is smaller than for freely diffusing colloids 
{a = 1/2), the dipole field of the COM of the subchains must decay more slowly and one 
expects, accordingly, a much weaker decay of the associated VCF. 



3. C enter- of-mass velocity correlation function 

Following the numerical strategy used in Sec. IIII Al for BFM beads, we characterize now 
in more detail the dynamical correlations already visible from the chain MSD h^{t) and the 
subchain MSD hs{t) by computing directly their second derivative with respect to time, i.e. 
the displacement correlation functions C^{t) ^ d^h^{t)/2 (Figs. O [8] and [9]) and Cs(t) ~ 

The displacement correlation function C^{t,5t) = {u{t + Iq) ■ uito)) ^^J for the dis- 
placement vector u{t) = r-^{t + 6t) — r-^{t) of the chain COM r-^{t) is shown in Fig. [7] for 
two chain lengths, = 16 (top) and = 8192 (bottom). Averages are again performed 
over all configuration pairs {to,t + to) possible in the set of 10^ configurations sampled for 
each St indicated. As in Sec. IIII Al we find that C^{t, 6t) ~ 6t^ for t ^ 6t. For clarity, only 
the data subset is indicated which is used to construct the cummulated final VCF C-^lt) 
(as shown below in Fig. [8]). The bold lines represent the predicted short-time exponent 
u = 5/4. We emphasize that this exponent can be observed for A^ = 8192 over nearly five 
orders of magnitude. For A^ = 16 we also indicate the exponent a; = 5/2 expected for times 
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FIG. 7: VCF Ct^{t,6t) for two chain lengths, iV = 16 (top) and = 8192 (bottom), obtained 
using L26-moves. We only indicate for each 6t the data window which is used to construct (by 
logarithmic averaging) the final VCF CT<s{t). The bold lines represent the exponent a; = 5/4 which 
is generally observed for times 10 t <^ Tn- The exponent a; = 5/2 (dashed line) indicated for 
= 16 is expected for t ^> Tn- 



t ^ Tn where the chains should behave as colloids according to Eq. (j3]). The magnitude 
of the signal decreases strongly with N, which together with the fact that fewer chains per 
box are available makes the determination of the VCF more difficult with increasing chain 
length. 

The A^-dependence of the VCF for t ^ Tn is further analysed in Fig. |8] where we plot 
the reduced VCF y = —Cf^(t){b^p)/(WDjsf) vs. the reduced time x = Wt using the effective 
monomer mobility W and the diffusion coefficient -Dn determined above. This scaling makes 
the axes dimensionless and rescales the vertical axis by a factor A^. As shown by the 
successful data collapse for chain lengths ranging from A^ = 16 up to A^ = 8192 on the 
power-law slope indicated by the bold line, the VCF scales exactly as CN{t) ~ 1/A^ for 
t ^ Tn. This confirms the already stated "locality" of the correlations [37|, Eq. (fT4|) . i.e. 
the effective forces acting on subchains add up independently to the forces acting on the 
entire chain. We have still to motivate the precise form used for the rescaling of the axes. 
According to Eq. ([7]) we claim that the VCF C^{t) scales as a function of the reduced 
time t/T-^. Substituting the typical chain size i?N ~ hN^/'^ and the chain relaxation time 
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FIG. 8: Reduced VCF y = -C]<i{t){b^p)/{WDj<f) for different chain lengths as a function of 
the reduced time x = Wt with W being the local mobility, ^ 1/N the self-diffusion constant. 
For times t ^ Tn the VCF thus scales as CN(t) ^ ^/N, i.e. the correlations must be due to ~ A'" 
local independent events. The exponent cj = 5/4 (bold line) is observed over up to five orders of 
magnitude in time. 



Ttss ^ R^/ Dt<{ ^ /W and using w = 5/4 this yields 



y = -b^p 



{Wty^/^N^ for t < Tn 



(16) 



with c being an empirical dimensionless constant. The bold slope indicated in the plot 
corresponds to a value c ~ 1. Interestingly, since CN(t) ~ dlh^{t)/2, it follows from Eq. flT6]) 
that 



h^(t) = 6D^t 1 



16c 



my 



-1/4 



(17) 



As may be seen from Fig. [6] (stars) for = 8192, Eq. ( |T7I) with c = 1 provides an excellent 
fit of the measured h-^{t). We also note that the second term in Eq. f|T7|) dominates the 
dynamics for t t* = I^~^(16c/96^p)^ A^*^ ~ 10'^. This is indicated by the vertical arrow 
in Fig. [6l Hence, for t <^t* the stars correspond to an exponent P = 2 — u = 3/A, Eq. ([9]), 
being close to the phenomenological exponent /3 ~ 0.8 which motivated our study. The 
central advantage of computing the VCF C-^{t) lies in the fact that it allows us thus to make 
manifest that (negative algebraic) deviations from the Rouse behavior exist for all times and 
not just for t < t*. 

Returning to our discussion of Fig. [8] we emphasize that the VCF of shorter chains decays 
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FIG. 9: Collapse of the reduced VCF y = C^{t) /Ct^ with Ct^ ~ 1 /N'^^'^ as a function of a; = t/T^ 
as suggested by Eq. ([7]). The large time behavior (x 1) is described by the exponent w = 5/2 
(dashed line) expected for overdamped colloids, Eq. ([3]). The exponent oj = 5/4 (bold line) for 
X ^ 1 is confirmed over nearly eight orders of magnitude. This is the central numerical result of 
this study. 

more rapidly for large times following roughly the dashed line corresponding to the exponent 
= 5/2 expected for effective colloids. We verify in Fig. [9] that this bending down of the 
VCF is described by the announced scaling in terms of a reduced time x = t/T^ and a vertical 

axis y = CN{t)/CT^ using the amplitude Ct^ = C^^it = Tn) ^ -{R^/T^fp*/p l/A^^/^ 

stated in Eq. ([7]). As shown by the data collapse, there is only one relevant time scale in 
this problem, namely the chain relaxation time Tn, for which the deviations for short times, 
where polymer physics matters (w = 5/4), and for large times, where polymer chains behave 
as effective colloids {uj = 5/2), nicely match. It is worthwhile to emphasize that Eq. ([7]) 
together with the locality of the deviations, C^{t) ~ 1/A^, immediately imply the exponent 
u. This can be seen by counting the powers of the chain length N, 

-1 = {1/2-2)2 + {l-d/2) + 2uj, (18) 

which implies u = ((i + 2)/4 = 5/4in agreement with Eq. and the numerically observed 
time dependence. Assuming Eq. ([7]), the exponents for and t thus contain the same 
information. 

The subchain VCF Cs{t) presented in Fig. [10] has been obtained as the total chain VCF 
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FIG. 10: VCF Cs{t) for subchains of arc-length s obtained for chains of length N = 1024 using an 
overlap penalty e = 10 and L26-moves. Main panel: For t ^ Tg ~ the subchains scale as the total 
chains, i.e. Cs{t) ~ s~^t~^/^ (bold line). For intermediate times Tg ^ t ^ Tn the displacement of 
subchains (1 ^ s ^ N) follows the monomer MSD and we find thus Cs(t) ~ sO(6VF)2(m)V2-2 
(dash-dotted line). Inset: Setting Tg = s'^/W and Ct, = -{bWfs-^ and plotting y = Cs{t)/C^^ 
as a function of x = t/Tg the data scales for x » 1. The dashed slope corresponds to y = l/8x'^/^. 



CN(t), the only difference being that the COM rs{t) of the subchain defines now the displace- 
ment vector u{t) = rs{t + 6t) — rs{t). The data presented in the main panel is rescaled as in 
Fig.[8]with Dg = 0.009/s setting now the diffusion constant. For small times t ^ Tg ~ s'^/W, 
all data sets collapse on the same slope with exponent u = 5/4 (bold line) as in Eq. (fTBIl . 
This confirms that 

in agreement with Eq. ( fl4l) . This shows that the same deviations occur for arbitrary sub- 
chains and that the colored forces acting on the subchains add up independently to the 
effective forces acting on the total chain. Since for intermediate times Tg ^ t ^ the sub- 
chains are enslaved by the monomer motion according to Eq. f|T5|) . it follows using Eq. ( ITT]) 
that C,{t) ~ s°{bW)\Wty/^-^ (dash-dotted line). This scaling can be better seen in the 
inset of Fig. [TD] where we plot y = Cs(t)/CT^ as a function of a; = t/Tg setting Tg = s'^/W 
and = —{bW)^s~^. This allows to bring to a nice collapse the subchain VCF for x ^ 1. 
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x = Wt 

FIG. 11: VCF CN(t) for various monomer overlap penalties e for chain length N = 1024, volume 
fraction </> = 0.5 and L26-moves. For sufficiently large times (t <C Tn) all data appears to collapse on 
the same asymptotic exponent = 5/4 (bold line) for incompressible polymer solutions. Additional 
physics is visible at short times. 

The dash-dotted hne corresponds to the expected decay y = l/8x^/^ for intermediate times. 
C. Robustness of scaling behavior 

We have focused above on one specific choice of operational parameters. In order to check 
whether our key predictions Eq. ([7]) and ([8]) remain valid more generally under different 
conditions, additional simulations have been performed which we summarize here: 

(i) Switching to the L06-moves of the original BFM algorithm, but keeping the same 
static conditions (0 = 0.5, e = 10), the same scaling is obtained for the MSDs and the VCFs 
as before (not shown). Using the directly measured local monomer mobility W = 0.002 and 
the self-diffusion coefficient = 0.007/A^ we get for the VCF exactly the same scaling as 
in Fig. IHl This is consistent with the fact that Eq. ([7]) does not depend explicitly on the 
specific local dynamics. 

(a) Using again L26- moves we have varied the overlap penalty e at constant chain length 

= 1024 and volume fraction (p = 0.5. Since Eq. ((Tj) does only depend implicitly on g, this 
suggest that Eq. (fT6|) for times t ^ Tn should also remain valid if one uses for rescaling 
of the axes the measured effective bond length b{e), the mobility W{e) and the diffusion 
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constant D-f<s{e) indicated in Tab. [H This is in fact borne out for the asymptotic behavior of 
the rescaled VCF displayed in Fig. [11] In analogy to the angular correlation function Pg for 
soft melts presented in panel (c) of Fig. [31 the deviations from Eq. fll6p for short times are 
expected since the incompressibility constraint is only felt if distances corresponding to the 
static screening length are probed. 

(in) A similar plot has been obtained for C-^(t) if we change the dimensionless compress- 
ibility by varying the volume fraction (p at constant chain length = 1024 and overlap 
penalty e = 10 (not shown). As long as the chains remain sufficiently entangled, all data 
merge for large times on the asymptotic behavior of incompressible melts, Eq. ( [TB|1 . for 
t < Tn. 



IV. CONCLUSION 



A. Summary 

The incompressibility constraint of polymer melts is known to restrict the fluctuations of 
sub) chains and generates thus scale-free static deviations from the Gaussian chain statistics 
12|-|l4|. In this paper we addressed the question of whether the incompressibility constraint 
also causes dynamical correlations of the (sub) chain displacements. To avoid additional 
correlations we have supposed that the dynamics is locally perfectly overdamped ("no mo- 
mentum conservation") and that the chains may freely intercept ("no reptation"), i.e. the 



relaxational dynamics is assumed to be to leading order of Rouse type [10| . We have shown 
that the above conditions may be realized computationally by a variant of the BFM using 
finite monomer excluded volume penalties and local topology non-conserving MC moves 
(Fig. [2]). Sampling chain lengths up to = 8192 allowed us to carefully check the A^-scaling 
of the deviations (Figs. [HI [9]). Such deviations are visible from the short-time scaling of the 
COM MSD h^it) presented in Fig. [6] confirming published computational and experimen- 



tal work 23N31|. We have shown that a better characterization of these deviations can be 
achieved by means of the COM displacement auto-correlation function Cf^(t) ^ d^h^(t)/2. 
Computing directly the curvature of hj<s{t) this removes the large free diffusion contribution 
to the chain motion and allows to focus on the correlated random forces which cause the 
second term contributing to h^(t) according to Eq. ( |T71) . How such a "velocity correlation 
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function" can be computed within a MC scheme has been first illustrated for dense BFM 
beads (Fig. [5]) confirming the negative algebraic decay expected for overdamped colloids, 
Eq. (I3D. The observed exponent uj = {d + 2)a is expected j2| due to the coupling of a tagged 
colloid to the gradient of the collective density dipole field (Fig. H]) decaying in time by free 
diffusion {a = 1/2). As shown in Fig. [HI the same exponents (dashed lines) characterize the 
motion of polymer chains for large times (t 3> Tn) where the chains behave as effective col- 
loids. More importantly, we have demonstrated that the observed short-time deviations for 
hf^(t), Eq. IQ, can be traced back to the negative analytic decay of the correlation function, 
CjJ^t) ~ -N-H^'^ for t < Tn with a; ^ 5/4 in ci = 3 (Fig. E]) in agreement with Eq. 
|34l |. That the correlation functions decays inversely with mass shows that the process is 
local js?!, i.e. the displacement correlations of subchains add up independently (Fig. [TU]) . 
Assuming according to the key scaling relation Eq. ([7]) the relaxation time Tn to be the only 
characteristic time, both asymptotic regimes can be brought to a data collapse (Fig. |9]). 
The short-time exponent w = (rf + 2)/4 is implied by Eq. ([7]) and the locality of the corre- 
lations, Eq. ( ITSl) . A deeper insight is obtained by generalizing the well-known displacement 
correlations of overdamped colloids (Fig. [1]) to the displacement field of subchains of length 
s ~ with 5t being the time window used to define the displacements. Since subchains 
repel each other due to the incompressibility constraint, a tagged subchain is pulled back 



to its original position by the subchain dipole field [38|. Since for times 5t ^ t ^ Tn the 
relevant dipole field decays much slower than for colloids (a = 1/2— 7-1/4), the correlations 
are much more pronounced. 

B. Outlook 

In this study we have deliberately tuned our model to avoid topological constraints. 
Obviously, these constraints are expected to matter for the dynamics of real polymer melts 



lOl, 



11| . It is thus of interest to see how the presented picture changes if topology is again 
switched on by using the topology conserving L06-moves of the classical BFM algorithm. 
Preliminary data comparing systems of chain length = 8192 is presented in Fig. [T2j The 
dynamics is still of Rouse-type at small overlap penalties as can be seen for e = 10. With 
increasing e the crossing of the chains gets more improbable and the topological constraints 
become more relevant. This can be seen for e = 30 (diamonds) and even more for e = 100 
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FIG. 12: Topology induced additional correlations in BFM melts for chain length = 8192. 
Main panel: MSD h^{t) vs. time t comparing L26- moves for e = 10 (spheres) to L06- moves for 
different overlap penalties e. Topological constraints become relevant for L06-moves with e S> 10. 
Inset: The reduced VCF —CT<}{t)N confirms that for all system classes the short-time dynamics is 
described by the same exponent w = 5/4 (bold line). 

(triangles). The latter data set clearly approaches the power-law slope 1/2 (dashed line) 
expected from reptation theory for times larger than the entanglement time Tg. The vertical 
arrow indicates a value for this time obtained from an analysis of h{t). Apart from the much 
larger chain length used, our data for e = 100 is consistent with the results obtained using the 



classical BFM with e = oo 



23 



26| . As in Fig. |6]the bold line represents the effective exponent 



/3 = 0.8, Eq. (|6]). Superficially, it does a better job for systems with conserved topology due 
to the broad crossover to the entangled regime. However, this apparent exponent is by no 
means deep as revealed by the displacement correlation function C-!ss(t) plotted in the inset. 
For short times all systems are well-described by the same exponent w = 5/4 (bold line) 
in agreement with Eq. ([8]). Note that the final points given for e = 100 decay slighly more 
rapidly. Unfortunately, neither the length of the analyzed time series nor the precision of 
our data does currently allow us to show that 

2 / , \ -3/2 



dc 



for t > Te 



(20) 
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Much longer 



1 /2 

as expected for reptating chains with ^ Te being the tube diameter 
time series are currently under production to clarify this issue. The numerical demonstration 
is obviously challenging, since the difference between the exponents 3/2 and 5/4 is rather 
small. In any case it is thus due to the dynamical correlations first seen in the BFM simula- 



tions of Paul et al. 



23 



3l| that the crossover between Rouse and reptation regimes becomes 



broader and more difficult to describe than suggested by the standard Rouse-reptation the- 
ory [l^, 11| which does not take into account the (static and dynamical) correlations of the 



composition fluctuations imposed by the incompressibility constraint 



39|. 
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